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a model  for  ductile  fracture  is  applied- €o  the  analysis  of 
failure  modes  on  a carbon-phenolic  composite,  loaded  at  the 
surface  by  a laser-supported  detonation  wave.  It  is  demon- 
strated that  for  a 100  joule /raj,  10  ns  pulse  the  model  pre 
diets  failure  of  the  sample  arthe  third  reflected  tensile 
stress  wave  propagates  through  it.  It  is  argued  that  high- 
intensity,  short  pulses  are  capable  of  producing  damage  * — ' 
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efficiently  in  composites,  but  that  with  LSD  surface  loading 
rather  than  direct  laser-target  interaction  the  rapid  stress 
relief  required  for  front-surface  spall  failure  is  difficult 
to  achieve.  The  attenuation  of  stress  waves  in  the  composite 
due  to  dispersion  and  viscous  dissipation,  the  effects  of 
defects  in  the  virgin  material,  and  the  nucleation  and  growth 
parameters  for  cracks  in  different  composite  materials  are 
critical  factors  in  the  analysis. 
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1.  INTRODUCTION 


The  general  objective  of  the  research  described  In  this 
report  Is  to  understand  how  composite  materials  may  be  damaged 
by  the  surface  stress  loading  resulting  from  the  Ignition  of 
laser-supported  detonation  (LSD)  waves  In  the  atmosphere  adja- 
cent to  the  surface.  This  type  of  surface  loading  Is  distin- 
guished from  others  (flyer  plate , x-ray , electron  beam,  direct 
laser-target  Interaction)  by  the  relatively  weak  mechanical 
coupling  between  the  region  in  which  the  stress  is  created  and 
the  bulk  of  the  target  material*  This  results  both  from  the 
mismatch  in  mechanical  impedance  between  the  air  and  the  solid 
target,  and  from  the  physical  separation  of  the  air  shock  and 
the  solid  surface.  As  a result  of  this  weak  coupling,  the 
peak  stress  developed  in  the  target  for  a given  energy  input 
is  comparatively  low.  Furthermore,  the  stress  relief  at  the 
surface  tends  to  be  relatively  slow,  but  a large  total  impulse 
may  be  transferred  to  the  target  over  an  extended  period  of 
time  at  very  low  stress  levels. ^ 

These  features  of  the  LSD-target  interaction  imply  that 
simple  empirical  damage  criteria  developed  to  correlate  ex- 
perimental results  obtained  with  other  types  of  loading  may 
not  be  useful  for  the  LSD  case.  We  have  therefore  employed  a 
more  complex  criterion,  essentially  that  developed  at  Stanford 
Research  Institute  (SRI)  for  ductile  fracture, ^ which  is 
based  on  detailed  models  of  the  nucleation  and  growth  of  de- 
fects, together  with  a large  number  of  microscopic  measurements. 


1.  Schriempf,  J.  T. , "Response  of  Materials  to  Laser  Radiation," 
NRL  Report  7728,  1974. 

2.  Seaman,  L. , T.  W.  Barbee,  Jr.  and  D.  R.  Curran,  "Dynamic 
Fracture  Criteria  of  Homogeneous  Materials,"  Air  Force 
Weapons  Laboratory  Report  A7WL-TR-71-156,  1971. 
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In  th«  absence  of  direct  results  on  LSD-induced  damage  in  com- 
posites, this  approach  appears  to  be  more  realistic  than  the 
simple  maximum-stress,  total-impulse  or  strain-energy  models 
employed  previously. 

The  general  features  of  the  LSD-composite  interaction 

outlined  above  also  have  implications  for  the  characteristics 

of  the  laser  beams  which  may  be  expected  to  induce  significant 

damage.  For  a given  fluence  incident  on  the  target,  say  100 
2 

joule s/cm,  one  increases  impulse  by  reducing  the  intensity 
(lengthening  the  pulse  time)  but  increases  peak  stress  by  in- 
creasing the  intensity  (shortening  the  pulse) . Both  the 
strain-energy  and  SRI  damage  criteria,  in  a general  way,  are 
more  dependent  on  peak  stress  than  on  impulse.  The  predicted 
deunage,  for  moderate  fluences  at  least,  has  a threshold  at 
some  minimum  intensity,  and  increases  with  intensity  above 
this  level. 

There  are  upper  limits  on  both  intensity  and  fluence 
for  laser  beam  propagation  through  the  atmosphere,  which 
depend  on  the  wavelength  of  the  laser,  the  number  and  size 
of  aerosol  particles  in  the  air,  the  time  profile  and  optical 
characteristics  of  the  laser  pulse,  and  the  propagation  dis- 
tance between  the  laser  and  the  target.  Here  it  will  merely 
be  remarked  that  intensity  alone  is  not  a proper  criterion 
for  propagation  of  pulsed  laser  beams  and  in  fact  for  pulses 
shorter  than  10  ns  the  limit  is  more  accurately  expressed  in 
terms  of  fluence.  For  sufficiently  low  fluences,  i.e.,  short 
pulses,  instantaneous  intensities  well  above  the  continuous- 

wave  (C-W)  breakdown  threshold  for  intensity  can  be  propagated; 

(3)  2 

one  measurement'  ' yielded  a threshold  of  800  joules/cm  for  a 

0.02  ns  neodymium  laser  pulse.  The  corresponding  intensity 

is  over  100  times  the  C-W  breakdown  threshold,  an  extreme 


I 3.  Ireland,  C.  L.  M. , A.  Yi,  J.  M.  Aaron,  and  C.  Grey  Morgan, 

"Focal-Length  Dependence  of  Air  Breakdown  by  a 20-psec 
Laser  Pulse,"  Appl.  Phys.  Letter  24,  175,  1974. 
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case.  Above  about  10  watts/cm  , of  course,  nonlinear-optics 
phenomena  may  occur  which  provide  a genuine  intensity  limit 
for  air  propagation,  but  surface  stresses  of  several  kilobars 
can  be  achieved  for  intensities  well  below  that. 

Multiple-pulse  laser  systems  provide,  in  principle,  the 
capability  of  delivering  high  total  fluences  to  the  target 
without  sacrificing  the  intensity  and  peak  stress  conditions 
needed  for  inducing  damage.  There  are  a number  of  problems 
in  the  multiple-pulse  systems  which  require  study,  but  we 
have  chosen  to  consider  first  the  mechanical  characteristics 
of  the  composite  target  which  provide  the  basis  for  any  evalua 
tion  of  laser  system  requirements.  In  particular,  answers  are 
needed  to  questions  such  as  the  following: 

1.  Under  what  conditions,  if  any,  can  front-surface 
tension  be  produced  by  LSD-wave  loading? 

2.  Under  what  conditions  can  damage  be  induced  by 
shear  loading,  as  for  example  at  the  edge  of 
the  irradiated  spot? 

3.  What  front-surface  loading  conditions  are  re- 
quired to  produce  a reflected  tension  wave  from 
the  back  surface  of  a target  capable  of  pro- 
ducing damage? 

4.  How  are  the  damage  characteristics  of  a com- 
posite related  to  the  composition,  structure, 
porosity,  etc.  of  the  material? 

5.  Assuming  that  the  damage  process  is  initiated, 
under  what  conditions  will  the  target  proceed 
to  complete  fracture? 

The  work  described  below  is  a first  step  only  toward 
the  resolution  of  some  of  these  questions  (in  particular  the 
last  three) . It  is  believed  that  further  progress  along 
these  lines  may  result  in  a much  improved  insight  into  the 
laser-composite  interaction,  and  provide  a valuable  guide  to 
future  experimental  and  theoretical  efforts. 


» 

Section  2 contains  a description  of  the  ductile  fracture 
I model.  Sections  3 and  4 discuss  the  material  equation  of  state 

and  the  artificial  viscosity  model  used  in  the  calculations/ 
which  in  turn  are  described  in  Section  5.  Some  conclusions  are 
discussed  in  Section  6. 
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2.  MODEL  FOR  DUCTILE  FRACTURE 


The  ductile  fracture  model  we  have  employed  is  a modified 

12) 

version  of  that  of  L.  Seaman , et.  al. ' ' In  the  SRI  ductile 
fracture  model,  nucleation  of  cracks  occur  at  a rate 

a - a 

• • no 

N » N exp  

o <j^ 

if  the  ambient  tensile  stress  a exceeds  a threshold  cr  (other- 
. no 

wise  N ■ 0) . Similarly,  cracks  already  present  grow  at  a rate 


if  a exceeds  a threshold  where  R is  the  crack  radius. 

go 

The  following  constants  are  given  for  two-dimensional  carbon 

(4) 

phenolic  composites. 

0.2  kbar 

2.1  x 10^°  cm3/sec 
5 x 10~3  cm 
0.38  kbar 
0 

2 

89.3  dyne-sec/cm 

In  our  implementation  of  this  model,  the  radius  R is  replaced 
by  the  fractional  void  volume  V ^ N R3.  The  ambient  stress  a 
is  related  to  V by  the  linear  relation 


4.  Shockey,  D.  A.,  L.  Seaman,  M.  Austin,  D.-  R.  Curran,  and  R.  A. 
Armi stead  (Stanford  Research  Institute) , "Repetitively  Pulsed 
Laser  Induced  Fracture  as  a Damage  Enhancement  Mechanism," 
Proceedings  of  the  1975  DOD  Laser  Effects/Hardening  Conference, 
NASA  TM  X-73,084,  Vol.  II,  p.  373. 
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a » ag  (1  - acV)  (1) 

(2) 

where  a„  is  a constant  taken  to  be  2.0  and  a , the  stress 
c s 

in  the  solid,  is  assumed  to  be  independent  of  V and  can  be 
evaluated  at  the  beginning  of  each  step  from  the  known  value 
of  a and  V.  A second  relation  between  c and  V is  the  dif- 
ferential equation 

m * ***  r £) + h - V v m 

• "3 

where  VQ  N0Ro  * Th*s  equation  states  that  the  total  change 
in  fractional  void  volume  results  from  the  production  of  new 
void  nuclei  and  from  the  growth  of  voids  already  present.  Since 
this  equation  is  linear  in  V,  its  analytic  solution  may  be  used 
together  with  (1)  in  an  iterative  root-finding  procedure  which 
determines  a time-advanced  value  of  V and  a corresponding  modi- 
fied value  of  o.  Since  the  terms  on  the  right  side  of  (2)  can- 
not be  negative,  the  procedure  can  only  increase  the  void  volume 
and  decrease  the  ambient  tensile  stress.  When  V for  any  mass 

element  exceeds  the  value  0.95/a  a fracture  flag  is  set.  In 

c 

this  calculation,  the  fractured  region  was  treated  as  essentially 
void,  able  to  reflect  both  tension  and  compression  waves. 


3.  EQUATIONS  OF  STATE 

The  equations  of  state  for  the  fiber  and  resin  constit- 
uents are  described  in  detail  in  Reference  5.  A brief  summary 
is  presented  below.  The  Mie-Gruneisen  form 

P * GpE  + (1  - j Gy)  PH  (y) 

is  assumed  for  p > pQ  (compressed  materials)  where  P is  the 
pressure,  G the  Gruneisen  coefficient,  E the  specific  internal 
energy , 


and  PR  (y)  is  the  Hugoniot  pressure.  For  carbon  fibers,  with 
pQ  * 1.59,  the  Hugoniot  is  represented  by: 

PH  = 73.6  y (kbar) , 0 < y < 0.3 

■ 71.58  - 403.62  y + 795.37  y2  (kbar), 

y > 0.3 


The  Gruneisen  coefficient  for  carbon  fibers  is  assumed  to  vary 
with  E.  With 

T - 0.24  E < 7.247  kJ/g 

- 0.667  E > 66.24  kJ/g 


5.  Gurtman,  G.  A.,  M.  H.  Rice,  R.  A.  Cecil  and  E.  W.  Sims, 
"Material  Response  Models  for  Radiation  Effects  on 
Advanced  Three-Dimensional  Composite  Heat  Shield  Mate- 
rials," Air  Force  Weapons  Laboratory  Report  AFWL-TR-74- 
144,  1975. 
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and  interpolated  linear  dependence  on  energy  for  intermediate 
values  of  E,  then 


<E>  - pi  / rdE 

o 


defines  the  Gruneisen  coefficient  for  carbon  fibers. 

For  phenolic  resin,  with  pQ  * 1.25  and  p > pQ, 


P„  * 69.6  y + 615  y2  - 839  y3  + 656  y4  (kbar) 

n 


G * G(p)  - ~ [0.52  + 0.68  y exp  (-1.238  y) ] . 


For  expanded  materials,  with  p < pQ,  let 


e ■ 1 » 


be  the  "engineer in? N strain  variable.  Then  the  pressure  is 
assumed  in  the  form 

po 

P « Gp  [E  - E (E)  (1  - exp  (N  e) ) ] 

9 P 

with 


o o s 


where  and  G.  are  the  bulk  modulus  and  Gruneisen  coefficients 
o o 

at  p - po. 


V M. 


i * 
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For  carbon  fibers  with  p < pQ, 

G - G(E)  [1  + e exp  (-  ^)  ] 

where  G(E)  is  defined  as  above  for  compressed  carbon  fibers. 
The  sublimation  energy  for  carbon  fibers  is  taken  as 

E (E)  - 7.247  kJ/g,  E < 7.247  kJ/g 

9 

- 58.99  kJ/g,  E > 66.24  kJ/g 

with  linear  dependence  on  E for  intermediate  values. 

For  phenolic  resin,  with  p < pQ 

G - G(p)  ■ 0.20  + 0.32  yp/p0 

and 

E (E)  » 836.8  (joule/g) 
s 

» J (E  - 836.8)  + 836.8 

- 2301.2 


E < 836.8  joule/g 
836.8  < E < 3033.4 
E > 3033.4 


For  composites  with  fibers  perpendicular  to  the  direction  of 
propagation,  the  shear  modulus  and  yield  strength  are  assumed 
zero.  The  fiber  and  resin  constituents  of  each  cell  are 
equilibrated  mechanically  and  thermally;  the  total  energy  and 
volume  of  a cell  is  divided  between  fiber  and  resin  so  as  to 
yield  the  same  pressure  and  temperature.  The  composite  is 
assumed  to  be  59.34  percent  fiber  and  40.66  percent  resin,  by 
volume,  for  the  perpendicular  fiber  orientation,  with  zero 
initial  porosity. 
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4.  ARTIFICIAL  VISCOSITY  MODEL 


Many  finite-difference  calculations  of  stress-waves 
employ  a quadratic  artificial  viscosity 

- -k^AulAu,  Au  < 0 

for  numerical  stabilization  of  strong  shocks*  plus  a linear 
term 


q2  - -k2pc  Au 

for  elastic  solids,  where  p is  the  density,  Au  the  velocity  in- 
crease across  a zone,  and  c the  sound  speed,  with  k^,  k2  being 
dimensionless  coefficients.  The  quadratic  term  q^  is  always 
switched  off  for  expanding  zones,  but  the  linear  term  may  be 
used  in  expansion,  compression  or  both.  The  linear  term  q2  is 
useful  for  stabilizing  the  numerical  calculation  for  weak  shocks 
in  solids,  since  the  q^  term  is  too  small  to  be  effective  for 
these  materials.  However,  the  q2  viscosity  is  dissipative  and 
may  attenuate  narrow  stress  pulses. 

A quadratic  formulation  of  artificial  viscosity  for 
elastic  solids was  therefore  introduced: 

q3  - -k3  <Ap)2/p 

* -k3  Au  |Au  • |£{  , Au  > 0 


6.  White,'  J.  W. , "A  New  Form  of  Artificial  Viscosity  for 
Elastic  Solids,"  J,  Comp.  Phys.  16,  119-126,  1974. 
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where  p is  the  mean  pressure  in  the  shock,  and  the  second  form 
may  be  obtained  by  using  the  second  and  third  Rankine-Hugoniot 
relations.  The  ratio  Ap/AE  was  evaluated  as 


dP 

WE 


/ 


dE 

ca- 


using single-step  time  differencing.  The  term  was  used  only 
for  expanding  zones.  In  the  calculation  described  later  in 
this  report,  it  was  successful  in  calculating  the  initial  com- 
pressive pulse  propagation  with  much  less  attenuation  than 
occurred  with  q2,  but  failed  to  allow  a stable  reflection  from 
the  back  surface  of  the  composite  layer.  The  calculation  was 
therefore  continued  from  this  stage  with  q^  and  q2  terms  only. 
The  instability  of  the  q3  model  at  the  back  surface  is  not  at 
this  point  fully  understood,  but  it  should  be  noted  that  the 
model  was  developed  for  homogeneous  materials  with  simple 
elasticity  rather  than  for  the  more  complex  material  in  our 
application. 

The  spreading  of  the  wave  profiles  resulting  from  the 
q2  formulation  is  too  severe  to  be  totally  ignored.  While  in 
the  present  series  of  calculations  the  energy  dissipation  due 
to  q2  is  numerical  in  nature,  it  should  be  noted  that  most 
polymeric  materials  (such  as  phenolic)  are  somewhat  viscous. 
This  time  dependence  is  not  currently  accounted  for  in  our 
rather  simple  solid  equation  of  state  for  phenolic  (see  Sec- 
tion 3).  The  calculated  pulse  spreading  and  resultant  attenua- 
tion might  therefore  be  a reasonable  approximation  to  2D CP 
constitutive  behavior,  but  it  would  be  necessary  to  compare 
numerical  results  with  experimentally  determined  wave  profiles 
to  ascertain  this  to  with  any  degree  of  confidence.  We  might 
note  also  that  the  geometric  (and  non-dissipative)  dispersion 
inherent  in  any  laminated  composite  made  up  of  constituents 
possessing  differing  shock  impedances,  is  correctly  modeled  in 
the  existing  code  formulation. 
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5.  CALCULATION  RESULTS 

The  prototype  calculation  here  presented  is  intended 
to  exhibit  the  performance  of  the  ductile  fracture  model  and 
the  other  analytic  techniques  which  have  been  developed  for 
the  study  of  laser-composite  interactions.  Additional  calcu- 
lations of  this  type  would  be  required  in  order  to  infer 
actual  failure  criteria  for  the  material  represented , since 
the  parameters  of  the  calculation  were  purposely  chosen  to  be 
relatively  favorable  to  the  production  of  damage. 

The  target  material , a two-dimensional  carbon-phenolic 

composite,  was  described  in  Section  3 above.  The  thickness 

of  the  sample  was  0.254  cm  (100  mils).  The  incident  laser 

pulse  fluence  was  100  joules/cm  and  the  total  pulse  time 

10  ns.  The  arbitrarily  chosen  pulse  shape  (shown  in  Figure 

1)  and  wavelength  (1.06  ym)  had  very  little  effect  on  the 

calculation;  the  neodymium  wavelength  was  used  since  a C02 

laser  pulse  with  this  fluence  and  length  would  not  propagate 

through  clean  air.  The  reflectivity  of  the  surface  was  unity, 

4 2 

and  the  density  of  LSD  wave  nucleation  sites  10  /cm  . The 
absorption  of  laser  radiation  in  the  LSD  wave  was  represented 
by  the  model  described  in  Reference  7.  One-dimensional  propa- 
gation was  assumed,  so  that  no  lateral  rarefactions  contributed 
to  the  stress  relaxation  at  the  surface  following  the  pulse. 

Figure  1 shows  the  peak  pressure  in  the  LSD  wave  for  the 
first  10  ns  and  the  ensuing  blast  wave  for  an  additional  15  ns. 
Figure  2 shows  the  pressure  in  the  air  adjacent  to  the  solid 
surface  for  the  same  initial  25  ns.  The  maximum  at  15  ns  re- 
presents the  arrival  of  the  blast  wave  resulting  from  the  cut- 
off of  the  laser  pulse.  The  surface  air  pressure  is  also  shown 
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Figure  1.  Peak  stress  in  air  shock  versus  time. 
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in  Figure  3 on  a longer  time  scale,  to  1.2  us.  This  very 
gradual  relaxation  effectively  prevented  any  tensile  stresses 
from  developing  near  the  front  surface.  In  a real  situation, 
lateral  rarefactions  and  radiative  cooling  would  provide  a 
more  rapid  relief  of  the  front- surface  stress,  at  least  at 
late  times. 

Figure  4 shows  the  wave  loci  and  damage  contours  on  a 
space-time  diagram,  with  compressive  waves  shown  as  broken 
lines.  The  first  tension  wave  produces  most  of  the  micro- 
cracks, but  is  incapable  of  driving  them  to  significant  size. 

The  contour  for  10”  3 void  volume  fraction  coincides  with  this 
wave  through  most  of  the  composite,  excluding  only  thin  layers 
near  the  front  and  back  surfaces.  The  second  tension  wave 
brings  the  void  fraction  to  order  10-1  throughout  the  central 
two-thirds  of  the  target.  The  curvature  of  the  wave  locus  re- 
sults from  the  attenuation  of  the  wave,  here  in  a nonlinear 
portion  of  the  Hugoniot. 

The  third  tension  wave  initiates  the  fracture  locus  at 
a point  about  75  percent  through  the  target,  at  time  4.25  us. 
From  this  point  on,  the  waves  reverberate  between  the  fracture 
boundary  and  the  surfaces  of  the  target. 

The  damage  contours  expand  under  tensile  stress  only, 
so  are  somewhat  irregular.  At  8 us  the  fractured  region 
occupies  most  of  the  interior  of  the  layer  and  is  still  ex- 
panding slowly. 

Figure  5 shows  the  peak  stress  in  the  first  and  second 
waves.  The  remaining  figures  are  plots  of  the  spatial  wave 
shapes.  Figure  6 shows  the  first  compression  wave.  Figures  7 
and  8 the  first  tension  wave.  Figures  9 and  10  show  the  second 
wave  and  Figure  11  the  third  compression  wave.  In  Figure  12 
the  third  tension  wave  has  just  been  interrupted  by  the  fracture 
of  the  material,  and  the  two  sections  have  reflected  as  separate 
compression  waves  in  Figure  13. 


COORDINATE  (g/cm 


U«t|  «i  Vi  MAii  CVClt  NU • • ijo  Tint*  &•*** 

li>  l «... 


10**l  • I S*Q  VS  HASS  CVCtC  NO..  *30  T I HE*  1.00*7-0* 


22 


Figure  7.  Stress  (dynes/cm  ) versus  mass  depth  (g/cm  ) at  1.00  ys. 
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Figure  8.  Stress  { dynes/cm  ) versus  mass  depth  (g/cm' 
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Figure  10.  Stress  (dynes/cm  ) versus  mass  depth  (g/cm 


Figure  11.  Stress  (dynes/cm  ) versus  mass  depth  (g/cnT)  at  3.81  ys 


Figure  12.  Stress  (dynes/cnr)  versus  mass  depth  (g/cnr)  at  4.41  ys 


6 . CONCLUSIONS 


The  laser  pulse  chosen  for  the  calculation  described 
above,  100  joules/cm  , in  10  seconds,  is  clearly  well  above 
the  minimum  required  for  fracture  of  the  target.  A lower 
fluence,  a longer  pulse  time,  or  a thicker  target  would  be 
required  for  an  investigation  of  thresholds. 

The  calculation  underlines  the  critical  importance  of 
void  nucleation  in  damage  processes  described  by  the  ductile 
fracture  model  we  have  used.  Since  void  growth  is  exponential 
and,  for  this  material,  has  a zero  stress  threshold,  it  might 
at  first  appear  that  almost  any  stress  loading  capable  of 
nucleating  cracks  would  automatically  produce  ultimate  fail- 
ure. A closer  inspection  of  the  behavior  shows,  however,  that 
this  is  not  true  at  lower  stresses.  Near  the  front  and  back 
surfaces,  where  wave  overlap  limits  the  magnitude  of  the  ten- 
sile stress  to  values  near  threshold,  the  cracks  grow  very 
slowly.  The  fracture  model  near  threshold  is  therefore  de- 
pendent upon  total  impulse  as  well  as  upon  peak  stress.  With 
repeated  pulses,  as  remarked  above,  the  required  impulse  could 
be  supplied  over  a period  of  time. 

The  calculation  did  not  assume  initial  porosity,  since 
the  empirical  parameters  used  in  the  fracture  model  do  not 
specify  it,  and  arbitrarily  inserting  a porosity  would  there- 
fore have  been  inconsistent.  A refinement  of  the  model  is 
needed  to  represent  initial  defect  distributions,  which  may 
strongly  influence  failure  behavior. 

Finally,  the  calculation  demonstrates  the  effectiveness 
of  high  laser  pulse  intensity,  in  contrast  to  high  fluence  in 
producing  damage,  and  provides  a basis  for  further  investiga- 
tions of  damage  phenomena. 
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